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ABSTRACT 


The  general  differential  equations  describing  the  behavior  of  the 
rotary  regenerator  including  longitudinal  heat  conduction  in  the  direction 
of  fluid  flow  are  sufficiently  complicated  to  preclude  a  complete  analyti¬ 
cal  solution.  A.  fev»  solutions  found  in  the  literature  restricted  to  very 
special  cases  are  discussed. 

A  numerical  finite  difference  method  is  presented  which  will 
determine  the  effect  of  longitudinal  heat  conduction  in  rotary  regenerators 
for  steady  state  conditions.  In  the  development  no  assumptions  are  made 
which  would  restrict  the  range  of  parameters  for  which  the  analysis 
would  be  applicable. 

The  conduction  effect  on  the  regenerator  effectiveness  was  evaluated 
with  the  view  of  obtaining  results  most  useful  for  the  gas  turbine 
regenerator  problem,  however,  these  results  may  also  be  used  for  other 
regenerator  problems. 

A  CDC  1604  digital  computer  was  used  to  carry  out  the  computations. 

The  results  are  presented  graphically  and  in  tabular  form,  employing  a 
suitable  set  of  non  dimensional  parameters.  The  range  of  parameters 
which  have  been  covered  arc: 

0-5  Ciain/pBiax  =^1 -0 
1.0  ,  Cr/Cn^lO.O 
1.0  NTUo  4  20.0 
0.25  i  (hA)*j  1.: 

0  \  .  0.2 

0.15  As*  1.3 

The  author  expresses  his  appreciation  to  C.  P.  Howard,  Associate 
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I.  Introduction 

The  regenerative  type  of  heat  exchanger  has  been  subjected  to 
extensive  treatment  within  the  last  thirty  years  or  so  and  contributions 
to  the  theory  of  regenerators  have  been  made  by  many  writers.  The  most 
questionable  of  the  idealizations  made  in  the  development  of  the  theory 
i3  that  of  zero  matrix  conductivity  in  the  direction  of  fluid  flow.  The 
loos  due  to  heat  conduction  in  the  matrix  material  is  specially  important 
in  the  case  of  high  performance  regenerators  in  which  a  small  reduction 
in  the  effectiveness  will  result  in  a  large  increase  of  the  heat  transfer 
areas  to  compensate  for  the  loss.  Another  important  factor  is  the  modern 
trend  toward  smaller  sizes,  since  the  loss  due  to  conduction  of  heat  in 
the  matrix  obviously  increases  with  decreasing  length. 

r  1 

Schultz  2  solved  the  problem  including  conduction  for  the 
special  case  of  the  balanced  regenerator  and  very  large  values  of  Cr/C^n 
(i.e.,  the  capacity  rate  ratio  of  rotor  matrix  to  minimum  fluid  capacity 
rate).  In  this  case  the  system  of  partial  differential  equations  defining 
the  problem  for  the  steady  state  are  reduced  to  ordinary  differential 
equations  and  the  problem  can  be  solved  analytically. 

Hahnemann  3  extended  the  solution  to  the  unbalanced  regenerator 
(i.e.,  unequal  partitions  of  the  cross  sections  and  unequal  heat  transfer 
coefficients  together  with  unequal  capacity  rate  ratios  of  the  fluid 
streams)  but  still  for  an  infinitely  large  value  of  Cr^^in* 

London  has  shown ,  as  reported  by  Lambertson  i  ,  that  a  correction 
factor  for  conductivity  in  the  direction  of  fluid  flow  for  approximately 

^■Number  in  brackets  refers  to  bibliography  on  page  2.1 
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equal  fluid  capacities,  of  the  form: 

/L\  tg  _  i-t  A 

C  L 

should  result  in  a  somewhat  pessimistic  prediction  of  the  reduction  in 
effectiveness. 

It  is  the  objective  of  this  thesis  to  develop  and  present  a  finite 
difference  numerical  analysis  with  solutions  to  solve  for  the  conduction 
effect  on  effectiveness  for  the  rotary  regenerator.  Since  there  arc  six 
dimensionless  parameters  required  to  specify  behavior  in  the  most  general 
case,  the  ranges  of  the  parameters  covered  in  the  solutions  will  be 
those  of  most  interest  to  the  gas  turbine  regenerator,  however,  the 
analysis  will  place  no  restrictions  on  these  parameters. 
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2.  Method 


The  conventional  idealizations  and  boundary  conditions  assumed  in 
the  derivation  of  the  governing  differential  equations,  including  the 
heat  conduction  in  the  direction  of  fluid  flow  4  ,  are  the  following: 

(See  Fig.  la) 

a.  The  thermal  conductivity  of  the  matrix  material  is  zero  in  the 
direction  of  matrix  metal  flow.  It  has  a  finite  value  in  the 
direction  of  fluid  flo\/,  and  is  infinite  in  the  other  direction 
normal  to  the  fluid  flow. 

b.  The  specific  heats  of  the  two  fluids  and  matrix  material  are 
constant  with  temperature. 

c.  Ho  leakage  of  the  fluids  occurs  either  due  to  direct  leakage  or 
carry  over,  and  each  fluid  flow  is  unmixed. 

d.  The  convective  conductance  between  the  fluid  and  the  matrix  are 
constant  with  flow  length. 

e.  The  fluids  pass  in  counter  flow  directions. 

f.  Entering  fluid  temperatures  are  uniform  over  the  flow  inlet  cross 
section  and  constant  with  time. 

g.  The  matrix  temperature  gradient  in  the  direction  of  fluid  flow  is 
zero  at  the  ends. 

h.  Regular  periodic  conditions  are  established  for  all  matrix 
elements,  i.e.,  steady  state  condition. 

The  differential  equations  derived  by  Schultz  \2  are  given  in 
Appendix  1. 

It  was  suggested  b\  IX’.sVa^rra  md  carried  out  by  Lambertson  1  for 
the  case  of  no  conduction,  that  a  rotary  regenerator  could  be  divided 
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element 


FIG.  la  ILLUSTRATIVE  MATRIX  ARRANGEMENT 


AND  FLUID  FLOW 


into  finite  elements  as  represented  schematically  in  Fig.  lb.  This  same 
method  can  be  extended  to  the  case  with  conduction,  if  the  elements  are 
considered  to  be  fixed  ir.  space,  each  element  can  be  regarded  as  a  cross 
flow  heat  exchanger  with  a  gas  stream  and  a  metal  stream.  For  an  element 
on  the  side  of  the  heat  transfer  rate  by  convection  is  equal  to  the 

rate  of  change  in  enthalpy  of  the  fluid  across  the  element, 


vl  ) 


N  V 


(1) 


where  _iTavg  represents  the  mean  temperature  difference  between  the  fluid 
and  the  matrix  element.  For  a  snail  enough  element  the  arithmetic  mean 
temperature  difference  may  be  assumed  valid,  so  that 
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where  the  fluid  and  matrix  temperatures  are  subscripted  to  indicate  the 
average  temperatures  across  the  inlet  and  outlet  of  the  element  for 
simple  unmixed  cross  flow. 

Considering  an  energy  balance  for  the  element,  the  energy  transferred 
to  the  element  by  convection  plus  energy  transferred  to  the  element  by 
conduction  less  the  energy  transferred  out  of  the  element  by  conduction 
must  equal  the  energy  stored  in  the  element: 
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FIG.  lb 


SCHEMATIC  REPRESENTATION  OF  A  ROTARY  REGENERATOR 


By  using  equations  (l)  and  (2)  the  outlet  temperatures  for  an 

element  on  the  side  of  may  be  solved  for  in  terms  of  the  remaining 

temperatures  and  the  result  expressed  in  dimensionless  form  (Appendix  1) 
Since  the  longitudinal  heat  conduction  is  zero  at  the  ends,  the 
energy  balance  equation  is  different  from  equation  (2)  for  the  elements 
of  the  first  and  last  row.  Three  different  expressions  for  the  matrix 
outlet  temperatures  are  obtained: 

First  row  elements 
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i) 


Middle  row  elements 


Lsst  row  elements 


(3) 


(A) 


By  solving  for  the  fluid  outlet  temperature  on  the  side  of  in  terms 

of  the  matrix  outlet  temperature  from  equation  (1),  only  one  expression 

is  obtained,  which  holds  for  all  the  elements  on  the  side  of  C  , 
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By  the  same  method  the  outlet  temperatures  of  the  elements  on  the  side 
of  may  be  obtained  using  equations  (7)  and  (8). 

First  row  elements 
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Middle  row  elements 
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Last  vow  elements 


(U) 


Solving  for  the  fluid  outlet  temperature  for  the  elements  on  the  side 
of  from  equation  (7)  gives 


(12) 


which  holds  for  ail  the  elements  on  the  side  of  Cpiin.  For  the  values 
of  the  constants  see  Appendix  1. 

In  addition  to  the.  four  dimensionless  groupings  for  no  conduction, 
Cmin/Cmax  =  capacity  rate  ratio  of  fluid  flow  streams 
Cr/C^in  **  capacity  rate  ratio  of  rotor  matrix  to  minimum 
fluid  capacity  rate 


•J|f 

(hA)  =  conductance  ratio 


over-all  number  of  transfer  units; 


NTUo  = 


two  additional  dimensionless  groupings  are  needed  in  the  coefficients 
of  equations  (3),  (4),  (5),  (6),  (9),  (10),  (11),  and  (12)  to  take 
conduction  into  account.  These  conduction  parameters  may  be  defined  as: 
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conduction  parameter  on  side  of  C 


max 


conduction  parameter  on  side  of  Cpj^. 
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Also 


since  the  thermal  conductivity  and  the  length  of  the  matrix  are  the  same 
for  both  sides,  and  where 


c 


_  — — —  ,  conduction  area  ratio. 

A  x 

Another  form  of  defining  the  conduction  parameter  is 


_  - -  ,  total  conduction  parameter 

V.  (  , 

In  this  case  As  is  the  total  solid  area  available  for  conduction,  that  is 


This  form  offers  the  advantage  that  for  a  given  value  of  the  conduction 

parameter,  the  resulting  effectiveness  will  not  be  affected  by  small 

* 

changes  in  As  . 

From  the  schematic  representation  shown  in  Fig.  lb,  it  should  be 
noted  that  the  left  edge  of  the  regenerator  is  physically  the  same  os 
the  right  edge  and,  therefore,  the  matrix  inlet  temperature  for  the 
elements  of  the  first  colutan  on  the  side  of  is  the  same  as  the 

outlet  temperature  of  the  corresponding  element  of  the  last  column  on 
the  side  of  Expressed  mathematically. 


where  i.  and  f  can  take  the  values  of  l,  2,  3, . b'r. 

This  is  referred  to  as  the  reversal  condition. 

From  equations  (3),  (4),  (5),  (9),  (10),  and  (11)  It  is  seen  that 
in  order  to  solve  for  the  outlet  matrix  temperature  for  a  particular 
element,  it  is  necessary  to  know  the  matrix  temperatures  of  the  next 
element.  An  estimate  of  these  values  can  be  made  by  determining  the 


temperature  distribution  for  the  ease  of  zero  conduction  in  the  direction 
of  fluid  flow,  since  the  constants  D3  and  Dg  in  equations  (3),  (4)  and 
(5)  and  the  constants  F3  and  Fg  in  equations  (9),  (10)  and  (11)  are  zero 
in  this  case. 

A  temperature  scale  can  be  used  for  which  the  fluid  entrance  tempera 
ture  is  zero  at  Che  side  of  Cnjn  and  unity  at  the  side  of  Ctaax. 


r‘(i,g)  "  ° 

for  g  =  1,2,3, . 

. Kn 

x(l,j)  -  1 

for  j  =  1,2,3, . 

For  calculation  purposes  a  temperature  distribution  is  assumed  on 
the  left  edge  and  the  problem  is  solved  for  the  no  conduction  case  in 
order  to  obtain  an  initial  estimate  of  matrix  temperatures  for  working 
the  problem  with  conduction.  This  also  provides  the  no  conduction 
effectiveness  necessary  for  comparison. 

The  outlet  temperatures  are  calculated  for  every  element  by  repetitiv 
use  of  equations  (3),  (A),  (5)  and  (&)  depending  on  the  Location  of  the 
element.  The  calculation  is  started  with  the  first  element,  first 
column  on  the  side  of  ClLax  and  then  working  down  the  column.  When  the 
first  column  of  elements  is  completed,  the  second  column  is  calculated 
and  so  on  until  the  outlet  temperatures  of  all  the  elements  on  the  side 
of  Cmax  are  determined.  It  was  found  that  convergence  could  be  enhanced 
if  two  passes  per  column  were  made  before  proceeding  to  the  next.  one. 

The  reason  for  this  is  that  the  equations  for  the  outlet  temperatures 
contain  the  matrix  temperatures  of  the  next  element  which  are  necessary 
to  estimate. 

At  the  seal  represented  by  the  double  line  in  Fig.  lb, 


Tx(i,Nx  +  1)  -  Tn ( f , l ) 


since  the  matrix  outlet  temperature  for  the  elements  of  the  last  column 
on  the  side  of  Cmay  is  physically  the  same  as  the  inlet  temperature  of 
the  corresponding  element  of  the  first  column  on  the  side  of  C^n. 

The  3arae  method  is  applied  to  the  side  of  as  to  the  side  of 

Ctsjax*  If  the  temperature  distribution  assumed  on  the  left  edge  was 
correct  it  would  then  be  duplicated  on  the  right  (i.e.,  the  reversal 
condition  is  fulfilled).  If,  however,  this  is  not  the  case,  the 
resulting  temperatures  on  the  right  side  .Tn(ffun  are  now  used  on 

the  left  Tx(i,l)  and  the  procedure  is  repeated,  After  each  pass  (i.e., 
the  complete  calculation  of  a  temperature  distribution)  an  energy 
balance  is  made;  and  before  the  solution  is  accepted  for  a  particular 
set  of  parameters,  the  heat  balance  error  together  with  the  reversal 
condition  has  to  be  fulfilled  to  the  specified  accuracy. 

The  effectiveness  of  a  heat  exchanger  is  defined  (Appendix  2)  as 


From  the  conditions  of  the  problem  the  above  expression  for  effectiveness 
reduces  to 


The  heat  balance  error  was  computed  from 


1  .  v  AV  1-  JO  £ 
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That  t’  5t>).4it  i ’.'ill  ••O'lVC.Ti'C  (I.-,, 


«.  i<a  r':wr°iy  l  conditio-.,  n.t 


and  this  heat  balance  *rror  converge  to  zero)  will  depend  on  the  number 
o £  elements  used.  More  Important  than  this  is  the  relation  existing 
between  the  numb.v  of  subdivisions  for  the  three  streams.  The  sufficient 
but  not  necessary  conditions  for  convergence  are  (Appendix  1)  that 


N\ 


\  A  *  ' 


\\ 


Y 


r. 


N 


Fig.  2  is  a  plot  of  the  different  values  of  effectivencoo  obtained 
for  the  following  set  of  parameters- 


C  .  /C  -  1.0,  Cr/CL,„  »  2.0,  (hA)*  -  1.0,  NTUo  -  10. 0  \=0.1,  As‘  =  1. 

Several  curves  are  obtained  depending  on  the  relation  existing  between 
the  number  of  subdivisions,  for  every  stream.  To  extrapolate  the 
effectiveness  to  an  infinite  number  of  elements,  these  curves  nay  be 

f 

approximated  by  straight  lines.  The  error  introduced  will  depend  on 
the  number  of  elements  used  and  will  not  be  constant  for  all  the  range 
of  the  parameters.  As  on  illustration  consider  the  sot  of  parameters 
used  to  plot  Fig.  2.  Three  cases  arc  investigated: 

1)  HJt  -  Kn  -  Nr/2 

2)  -V.  '  Nn  «  Mr 

3)  -!x  -  ?IU  -  2Hr 

A  linear  extrapolation  Is  performed  using  the  values  of  effectiveness 

obtained  for  values  of  1  -"  •  From  the  results  it  is  found 

Nr(Nx  +  Nn) 

that  the  minimum  Hr  necessary  to  give  an  error  in  the  cxtropolatcd 
effectiveness  of  less  than  0.01%  is  18,  for  which  =* 


76.57% 


(\J 

to 

•H 


1 A 


l/Nr(Nx+Nn) 


for  the  three  cns<  s.  The  determining  faetor  in  the  extrapolation  is 
the  nunber  of  matrix  elements.  The  higher  the  value  of  '.r  U3;  d  for  the 
extrapolation  the  smaller  the  error  in  the  straight  line  approximation. 
Unfortunately,  an  increase  in  Mr  will  bring  a  higher  increase  of  Kx  and 
Nn,  since  the  convergence  condition  for  this  example  is 

N 

\»  * 

That  this  is  only  a  sufficient  condition  can  be  seen  from  the  case  of 

Nx  -  =  Nr/2 

for  which  the  solution  converges  for  values  of  Nr  >5;  and  when 

Hx  -  Hn  -  Nr 

the  solution  will  still  converge  for  Nr>20.  It  is  true,  however,  that 
the  closer  the  number  of  matrix  elements  are  to  the  sufficient  convergence 
condition,  the  fewer  passes  it  takes  for  the  solution  to  converge  to 
the  specified  error  in  the  effectiveness.  Table  1  shove  the  number 
of  passes  necessary  to  obtain  a  .  olution  with  a  heat  balance  error  of 
less  than  0.0057,. 

For  a  given  set  of  parameters  the  effectiveness  vas  determined 
first  using  one  set  of  subdivisions  and  then  recalculated  increasing 
the  number  of  subdivisions.  A  linear  extrapolation  was  performed  to  obtain 
the  effectiveness  corresponding  to  an  infinite  number  of  elements.  The 
accuracy  for  each  value  of  effectiveness  was  specified  to  be  of  four 
significant  figures.  When  these  values  were  extrapolated  to  an  infinite 
number  of  elements,  the  accuracy  vas  reduced  to  only  three  significant 
figures  since  a  compromise  had  to  be  made  with  the  time  available  for 
computations  and  th<.:  number  of  runs  to  be  made.  However,  the  method 


TAolA  l 

Number  of  passes  necessary  to  obtain  an  effectiveness  with  a  heat 
balance  error  of  less  than  0.005*  for  the  following  at  t  of  parameters : 
<«ln/Cpax  -1.0  .  Cj/Cn^,  -  2.0  ,  <hA)*  =  1.0  , 


HTUo  -  10.0  ,  X*  °‘l  »  As*  -1.0 


Hr 

V 

1  X 

*»n 

Mo.  of  passes 

vV 

7, 

12 

12 

12 

27 

76.581 

12 

24 

24 

22 

76.581 

14 

7 

7 

41 

76.578 

14 

14 

14 

34 

76.576 

14 

28 

28 

27 

76.576 

16 

16 

16 

41 

76.574 

16 

32 

32 

32 

76.574 

18 

9 

9 

62 

76.573 

1.8 

18 

13 

50 

76.572 

10 

36 

36 

38 

76.572 

20 

20 

20 

60 

76.572 

j  22 

22 

22 

71 

76.571 

*  Ojo  is  obtained  by  a  linear  extrapolation  of  the  effectiveness 
for  Nc  and  Nr-1. 
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imposes  no  restriction  as  to  the  ultimate  accuracy  of  the  solution. 

When  working  with  very  high  values  of  Cr/C,n^n  the  computation  time 
becomes  considerably  long.  This  is  due  to  the  fact  that  the  final 
temperatures  differ  very  little  from  the  assumed  initial  temperatures  and 
a  large  number  of  elements  is  required.  The  starting  estimates  of  matrix 
temperatures  then  become  very  important.  For  the  limiting  case  of  Cr/Cjj,in 
equal  to  infinite,  there  is  no  temperature  variation  in  the  direction 
of  metal  flow  so  that  the  method  without  modifications  is  not  applicable. 
For  this  case,  however ,  the  exact  solution  of  Hahnemann  3  is  available. 

The  Fortran  system  was  used  to  program  the  above  method  for  the 
CDC  1604  digital  computer.  The  matrix  notation  employed  in  the  equations 
derived  previously  is  very  convenient  in  writing  up  the  program,  since 
several  loops  can  be  made  by  just  changing  the  subscripts  of  the  variables 
involved.  The  complete  computer  program  is  given  in  Appendix  2. 


1? 


3.  Discussion  o£  results 


The  solutions  generated  ace  presented  in  tabular  form  in  Tables  2 
through  25  ancl  graphically  in  Figs.  3  and  4.  The  accuracy  in  the 
effectiveness  was  specified  to  be  of  three  significant  figures.  It  was 
found  that  for 

0.7  (hA)*^  As*  ^1.3  (hA)* 

the  resulting  values  of  effectiveness  were  almost  identical  to  the  ones 
calculated  for 

As*  =  (hA)* 

so  that  these  were  the  values  used  in  preparing  the  Tables. 

By  using  N’fUo  and  X  os  defined  in  the  nomenclature,  which  contain 

(hA)  and  As'  explicitly,  the  influence  of  these  two  parameters  on  the 

•  -A  C? 

conduction  effect,  ,  were  found  to  be  very  small  for  values  of 
Cjj^n/Ctuax  ">  0.9  and  \^0.1,  which  is  the  range  for  the  gas  turbine 
regenerator  problem. 

Figs.  3  and  4  show  the  effect  of  conduction  on  effectiveness  for 
values  of  Cmin/Cmax  =  1.0  and  0.9  respectively,  with  different  values  of 
NTUo,  and  \  ,  up  to  a  limit  of  10%.  The  conduction  effect 

continues  to  increase  with  increasing  X •  For  small  values  of  A  this 
effect  increases  nearly  exponentially.  By  further  increasing  X  »  a 
saturation  is  reached  which  approaches  the  limiting  condition  of  \  =  oo. 

Fig.  5  is  a  representative  curve  whLch  shows  the  limiting  influence  on 
the  conduction  effect  by  continuously  increasing  X  for  a  given  set  of 
parameters . 

The  conduction  effect  for  values  of  Cr/C^n  =  10. 0  arc  almost  identical 
with  the  values  reported  by  Schultz  2  and  Hahnemann  [3  ]  ,  6  for  the 
limiting  case  of  =  oo. 

io 


l'or  email  values  of  NTUo  the  effect  of  increasing  Cr/Cciin  ^  to 
increase  slightly  the  conduction  effect,  as  can  be  seen  from  Fig.  6. 

For  higher  values  of  NTUo  the  effect  of  conduction  is  reduced  with 
increasing  Cr/C,^. 

The  effect  of  increasing  NTUo  upon  is  intimately  related  to  the 
value  of  Cp/Cnin  and  Cin£n/Craax.  This  effect  can  best  be  explained  by 
observing  Fig.  6.  For  the  case  where  C^n/C^ax  =1.0  the  effect  of 
increasing  NTUo  will  result  in  an  increase  in  4r  regardless  of  the  value 
of  Cr/C,n^n;  of  course,  the  higher  the  value  of  Cr/Cm^n  the  smaller  this 
effect  will  be  for  relatively  high  values  of  NTUo.  Now,  by  reducing 
^min/Snax  sn<*  increasing  Cr/G[|t^n  a  point  is  reached  in  which  an  increase 
in  the  NTUo  will  result  in  a  reduction  in  ^ .  The  point  where  this 
occurs  will  depend  upon  the  value  of  the  conduction  parameter  X  . 
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A. 


Surrw *& r y  * : ad  Co ncl  us  ions 

A  finite  difference  numerical  analysis  is  presented  which  v/ill 
determine  the  effect  of  the  longitudinal  heat  conduction  in  rotary 
regenerators  for  steady  state  conditions. 

The  ranges  of  governing  parameters  which  have  been  covered  are; 

0.5  ^^roin^max 
l.O^/C^JO.O 
1.0  6 iJTUo  ^20.0 
0.25  4(hA)*4l.O 
0  4  A  =,  0.2 
0.15  <As*^1.3 

The  ranges  of  parameters  of  interest  for  gas  turbine  regenerator 
application  are  plotted  in  Figs.  3  and  4.  A  more  accurate  determination 
of  the  conduction  effect  than  that  given  in  the  figures  can  be  made  by 
interpolating  values  from  the  Tables. 

The  computer  program  used  to  carry  out  the  calculations  on  the 
GDC  1604  digital  computer  is  included  in  Appendix  2.  This  program 
imposes  no  restrictions  on  the  parameters  for  v?hich  the  analysts  would  be 
applicable.  It  is  recommended,  however,  that  for  values  of  Cr/Cm^n>10 
the  solution  for  Cr/C„V£n  -  10  be  used.  For  this  case  the  computer  time 
for  a  solution  increases  considerably  and  the  error  introduced  by  assuming 
Or /Onin  =  10  is  negligible. 
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Conduction  effect 

> 

ro 

\ 

1 

> 

_ ) 

HTUo 

\ 

— 

0 

0.C15 

0.03 

~\06 

0.12 

0.24 

0.68  . 

j 

1 

1 

46.62 

.62 

1.17 

2.12 

3.64 

5 . 72 

3.05 

2 

60.08 

.85 

1.64 

3.04 

5.38 

P  £2 

13.03  j 

3 

66.75 

1.01 

1.95 

3.63 

6.46 

10.67 

i 

15.90  j 

4 

70. SO 

1.14 

2.19 

4.09 

7.26 

11.98 

t  T  C*"*  ■ 

n’.lj  j 

5 

73.79 

1.25 

2.40 

4.47 

7.91 

12.98 

1 

19.26  ; 

I 

1 

10 

81.19 

1.66 

3.15 

5.77 

9.98 

15.98 

23.2 

20 

36.62 

2.19 

4.07 

7.22 

12,05 

18.65 

26 . 35 

40 

90.46 

2.69 

4.92 

8.51 

13.77 

20.73 

s 

28.69 

100 

93.61 

3.24 

5.84 

9.85 

15.44 

22 . 60 

30 . 64 

500 

95.57 

3.51 

6.36 

10.66 

1 6 . 42 

23.76 

31 . 91 

1000 

95.83 

3.55 

6.43 

10.77 

16.62 

23.92 

32.10 
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U  ??  .  8  L040  f  TUDTT'M.  HIE.iT 

’'Vain  Await  =  1.0  ,  Cj  /‘'■min  ~  • 

('Ovmuc  m\; 

.o  ,  :w.)t 

;n  or 

*■  ~ 

.5 

1 

j 

e  % 

Conduction 

effect 

> 

% 

.  ) 

HTUo 

\ 

0 

.015 

.03  .06 

.12 

.21 

.24 

1 

49,11 

1.22 

3.79 

5.51 

6.75 

2 

64.92 

1.70 

5.62 

8.47 

10.64 

3 

72.81 

1.97 

0.  L  4 

10.10 

12.78 

4 

77.61 

2.16 

7.32 

H.17 

16.17 

5 

80.86 

2.30 

7.32 

11.93 

15.16 

6 

83.23 

2.42 

8.22 

12.55 

15.92 

7 

85.04 

2.51 

8 . 54 

13.03 

16.  52 

8 

86.47 

2.60 

8 . 81 

13.42 

17.01 

9 

87.63 

2.67 

9.04 

13.76 

17.42 

10 

88.59 

2.73 

9 . 24 

14.05 

17.76 

12 

90.11 

2.84 

9.57 

14 . 52 

1.3.32 

14 

91.25 

2.94 

9.83 

14.88 

18.76 

16 

92.14 

3.01 

10.05 

15.13 

19.11 

18 

92.86 

3.08 

10.23 

15.43 

1 9.40 

20 

i 

93.45 

3.14 

10 . 38 

15.63 

19.64 

500 

99.60 

2.12 

3.95  7.15 

12.29 

19.63 

28.54 

1000 

99.79 

2.17 

4.0  7.21 

12.37 

19.73 

28 . 67 
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TABUS  9  LONG ITUD1 HAL 

HEAT  CONDUCTION  EFFECT 

^In^Snox  “  l-° 

»  Cf/Prain 

-  10.0  , 

(hA)*  *  Aa* 

-  .5 

e  % 

Conduction  ef£ect; ~  % 

NTUo 

0 

.015 

.03 

.06 

.12 

.24 

.48 

1 

49.96 

0.65 

1.23 

2.24 

3.84 

6.04 

8.51 

2 

66.59 

0.90 

1.72 

3.20 

5.67 

9.32 

13.84 

3 

74.91 

1.02 

1.97 

3.70 

6.62 

11.04 

16.64 

4 

79.89 

1.10 

2.12 

4.00 

7.21 

12.09 

18.34 

5 

83.22 

1.15 

2.23 

4.21 

7.61 

12.80 

19.48 

10 

90.78 

1.28 

s 

1 

2.49 

4.73 

8.59 

14.51 

22.13 

20 

95.11 

1.36 

2.65 

5.05 

9.20 

15.58 

23.76 

i 
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TABLE  10  LOMGTTtTDISAT,  HEAT  CONDUCTION  EFFECT 

‘~min/Qnax  ~  1-0  ,  CjVCmin  ~  1.0  ,  (hA)  =  As*  -  .25 

e  % 

Conduction  effect  — 

.8 

NTUo  '■  A 

0 

.025 

.05 

.1 

.2 

.4 

1 

46.59 

1.02 

1.90 

3.38 

5.57 

8.25 

10.84 

2 

60.13 

1.43 

2.69 

4.88 

8.22 

12.50 

16.84 

3 

66.64 

1.70 

3.21 

5.80 

9.78 

14.88 

20.07 

4 

70.10 

1.92 

3.61 

6.50 

10.88 

16.47 

22.16 

5 

73.90 

2.10 

3.94 

7.05 

11.73 

17.65 

23.66 

10 

81.23 

2.74 

5.05 

8.81 

14.25 

20.93 

27.63 

20 

86.64 

3.52 

6.27 

10.56 

16.50 

23.61 

30.62 
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TABL 

E  11  LONGITUDINAL 

!!EAT  CONDUCTION  EFFECT 

! 

Saln^tax  ~  » 

Cf/Croin 

-  10.0  , 

(hA)*  =  As* 

=  .25 

1 

c  %  , 

1 

Conduction  effect  -^L  % 

NTUo  ^ 

0 

.025 

.05 

.1 

.2 

.4 

.8 

1 

49.96 

1.07 

1.99 

7.92 

5.84 

8.68 

11.49 

2 

! 

66.  60 

1.47 

2.79 

5.06 

8.59 

13.21 

18.05 

3 

! 

74.91 

1 

1.68 

3.20 

5. 84 

10.00 

15.53 

21.43 

4 

79.89 

1.80 

3.44 

6.32 

10.86 

16.93 

23.48 

5 

83.22 

1.88 

3.61 

6. 64 

11.44 

17.88 

24.86 

10 

90.78 

2.1 

4.03 

7.44 

12.83 

20.10 

28.02 

20 

95.11 

2.23 

4.29 

7.92 

13.69 

21.46 

29.83 
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7*.,  \2  LOr-.:*TUDI  :AL  \r6t.?  C  'T^-UCT iOl’  EFFECT 


\  \ 

J5Tt!u\ 


.  Vnin  *  ’  *  ’J  ,  V*> '  -  '«*  -  t  .0 

'*  % 


0 


Conduction  effect  -  —  c X 


,01 


.02 


.04 


.00 


.16 


l 

47.63 

0.38 

0.73 

1.36 

2.41 

3.98 

5.95 

2 

61.56 

0.53 

1.02 

1.93 

3.53 

6.  no 

9 . 62 

3 

68.46 

0.62 

I 

1.21 

2.30 

4 . 24 

7 . 37 

1 1 . 60 

*»4 

72.74 

0.70 

1.36 

2.60 

4 . 78 

8.32 

15.33 

5 

75.72 

0.78 

1.49 

7.84 

5.22 

9.06 

14.50 

10 

83.20 

1.01 

1.96 

3.70 

6.70 

11.39 

;  ? .  ?:> 

20 

88.69 

! 

1.32 

! 

i 

2.52 

4 . 66 

3.21 

13.56 

o  .*>  • 

♦  * 

40 

1 

92.41 

i 

1.68 

3.15 

5. 64 

9.58 

15.32 

22.74 

ICO 

95.50 

2.27 

4.08 

6.94 

11.17 

17.10 

24 .  67 
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Tn  HI. 

1 !.j: 

toiiynitft/ 1 

Hr,’AT  O' 

EFFECT 

**  * c* 

»  Cr/Gmin 

*2.0  , 

(hA)  *  =  A»* 

*  L.O 

-  % 

Conduction  effect^; —  % 

K‘ifUo\ 

0 

.01 

.02 

.04 

.03 

.16 

.32 

1 

i 

O 

o 

LI 

1 _ 

0.60 

0.77 

1.43 

2.53 

4.17 

6.23 

!  2 

66.92 

0.55 

1.07 

2.03 

3.71 

6.61 

10.16 

1  3 

75.30 

0.63 

L  .24 

2.37 

4.37 

7.66 

12.36 

\ 

80.40 

0.69 

1.35 

2.59 

4.82 

8.50 

13.80 

5 

83.85 

0.73 

1.43 

2.76 

5.14 

9.10 

14.82 

10 

) 

91.99 

0.85 

1.67 

3.23 

6.04 

10.71 

17.44 

20 

96.82 

0.90 

1.79 

3.49 

6.59 

11.74 

19.11 
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fABLi  14  LONGITUDINAL'.  HEAT  CONDUCTION  EFFECT 


^tin/Criax  3  *9  , 

^r^nin  ~ 

10.0  , 

(hA)*  *  As 

=  1.0 

*  % 

Conduction  effect  - —  % 

\ 

HTUo 

0 

.01 

.02 

.04 

.08 

.16 

.32 

l 

51.22 

0.41 

0.78 

1.45 

2.56 

4.23 

6.32 

2 

68.81 

0.56 

1.09 

2.06 

3.77 

6.49 

10.28 

3 

77.67 

0.63 

1.24 

2.38 

4.39 

7.68 

12.38 

4 

82.98 

0.68 

1.33 

2.56 

4.76 

8.39 

13.65 

5 

86.52 

0.70 

1.39 

2.68 

5.00 

8.87 

14.49 

10 

94.36 

0.75 

1.48 

2.88 

5.48 

9.87 

16.35 

20 

98.36 

0.65 

1.33 

2.70 

5.35 

10.03 

17.08 
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TABLE  15  LONGITUDINAL  !IEAT  CONDUCTION  EFFECT 


Cmin/(Wx  -  -9  , 

Cr/^rnin  = 

=  1.0  ,  (hA)* 

r  A  * 

==  As  = 

■  .5 

e 

r  o 

Conduction  effect 

)  '< 

» 

\ 

NTUo 

0 

.015 

.03 

.06 

.12 

.24 

.48 

1 

47.53 

0.54 

1.03 

1.88 

3.23 

5.11 

7.23 

2 

61.40 

0.74 

1.43 

2.67 

4.76 

7.89 

11.81 

3 

68.26 

0.87 

1.69 

3.17 

5.70 

9.55 

14.48 

4 

72.51 

0.97 

1.88 

3.55 

6.40 

10,73 

16.29 

5 

75.47 

1.06 

2.05 

3.86 

C  .94 

11.63 

17.64 

10 

82.98 

1.35 

2.61 

4.89 

8.69 

14.33 

21.38 

20 

88.34 

1.68 

3.23 

5.97 

10.39 

16.70 

24.35 
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tA  r  r 


If-  LONG 3 


L*tJDt;!AL  KEA* 


CONDUCT!  Oil 


V  Y  1 T '  V 

i'J-  r  '.t.  i 


r  /.•' 

"min  max  • 

9  »  T 

2.0  ,  (hA)*  = 

As*  =  .5 

e  % 

Conduction  effect 

)  € 

-  % 

s! 

NT  Go 

0 

.03 

.12 

,21 

.30 

1 

50.27 

1 . 03 

3.40 

4.96 

3.09 

66.83 

1.53 

5.08 

7.71 

9.73 

3 

75.23 

1.77 

6.01 

9.23 

11.76 

4 

80.32 

1.93 

6.63 

10.24 

13.08 

*1 

83.77 

2.05 

7.09 

10.97 

14.02 

6 

86.28 

2 . 14 

7.45 

11.53 

14.75 

7 

88.18 

| 

2.22 

7.74 

11.98 

15.32 

3 

89.69 

2.29 

7.97 

12.34 

15.78 

9 

90.90 

2.34 

8.17 

12.64 

16.17 

10 

91.91 

2.38 

8.34 

12.90 

16.49 

12 

93.47 

2.45 

8.60 

13.31 

17.00 

14 

94.62 

2.50 

3.80 

13.61 

17.39 

16 

95.50 

i 

j 

2.54 

8.96 

13.85 

17.67 

18 

1 

96.20 

2.56 

9.07 

14.03 

17.90 

20 

96.76 

2.57 

9.16 

14.18 

18.09 
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?Ar,U'  17  1.0' MI 

»  v  J/  i  .U  .L« 

UAT  COUfi'tCI 

ON  tFt 

aCT 

S'  !  - 

■*su«  T-- 

!<V'ln  51 

13.0  ,  CiA}* 

-  trr 

=  .5 

<2  % 

Conduction  effect 

) 

7o 

v 

\ 

OTUo\ 

\ 

0 

0.015 

0.03 

0.06 

0.12 

0.24 

1.48 

1 

51.  22 

9.58 

1.10 

2.01 

3.45 

5.45 

7.73 

2 

68 .  80 

0.80 

l .  35 

7 . 39 

5.15 

8.55 

12.83 

3 

77.66 

0.92 

1.78 

3.35 

6.04 

10.19 

15.54 

4 

82.98 

0.98 

1.91 

3.62 

6.58 

11.18 

17.22 

5 

86.51 

1.02 

1.99 

3.79 

6.39 

11.85 

18.33 

10 

94.36 

1.07 

2.12 

4.11 

7.65 

13.29 

20.77 

20 

93.35 

0.96 

1.96 

3.97 

7.71 

13.80 

21.91 
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T'VI, 

Silr/C»na 

o  13  1 O^C  riODINAL 

jt  =  •-*  »  1  r  /'in lu  ~ 

HtAT  CO».T*UCTIO:i  E? 

l.O,  (hA)*  =  As 

F£CT 

=  .25 

i 

1 

i 

j 

1 

e  % 

Conduction  effect 

,  ? 

_ 

OTUo 

0 

0.025 

0.05 

0.1 

0.2 

0.^ 

0.3 

1 

47.40 

O.G? 

1 . 62 

2.90 

4.62 

7.21 

9.57 

2 

i 

61.29 

1.18 

2.26 

4.13 

7.08 

10.98 

13.09 

3 

68.13 

1 .39 

2.65 

4.87 

t  .39 

13.10 

13.12 

4 

72.37 

1.54 

2.94 

5.41 

9.32 

14.54 

20.11 

5 

75.31 

1.66 

3.17 

5.83 

10.03 

15.61 

21 . 54 

10 

82.74 

1 

2.05 

3.92 

7.17 

12.16 

18.61 

25.35 

20 

i 

i 

33.02 

1 

1 - 

2.43 

4.67 

G.45 

14.04 

21.04 

20.16  ‘ 
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r  * 

,h*  ?  ^-i’C 

i  Lrl>  L  'hi L 

■  ffcl/iT  M  >)U 

TTiOt:  r.ri  i 

Cl' 

"'KiUv  1-4,..  *  > 

<\.  /< 

~  10  3  ,  ■ 

-*>  ^  T  As* 

=  .25 

e  % 

Conduction  1 1  f v c 

t,4~  % 

0 

f 

1  •i  ' 

>«  .  .  — 

0.05 

0.1 

0.2 

0.4 

n.s 

t’  ■*  > 

0.9  3 

1  4 

1.13 

5.  1  > 

*7  7  ' 

,  *  1  -f 

6' ,  no 

1 .29 

2,45 

4.49 

7. 72 

12. Ci 

16.74 

77 .  »'f 

1.47 

2.82 

".21 

°.06 

14.32 

20.  '2 

82.93 

1.58 

3.03 

r .  64 

>.37 

15.7C 

22.19 

00 

j-* 

1.64 

3.17 

5.93 

10.41 

16.62 

23.57 

9^.3'* 

1.7  4 

3.41 

6 „ 48 

11 .56 

18.67 

26.62 

93.33 

1.59 

3.26 

6.*  C 

1  < .  r  9 

19.54 

28.  rr 
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TA3LL  20  T.O: 

C’lUDF-VVL 

HEAT  ro 

■3>ucTio::  £r? 

EOT 

“  0.5 

>  Cr/Qntn 

-  1.0  , 

(IsA)  '  »  As* 

~  1.0 

e  %> 

Conduction  effect,-^—  % 

trruo 

0 

0.01 

0.02 

0.04 

0.00 

0.16 

2 

66.89 

0.35 

0.63 

1.29 

2.37 

4.15 

3 

74.01 

0.41 

0.80 

1.53 

2.84 

5.02 

4 

78.16 

0.45 

0.89 

1.71 

3.19 

5.65 

5 

80.89 

0.49 

0.97 

1.86 

3.46 

6.12 

10 

87.15 

0.65 

1  .26 

2.39 

4.38 

7.57 

20 

91.07 

0.86 

1.64 

3.04 

5.37 

8.94 

0.32 

6.71 

8.25 

9.29 

10.06 

12.15 

13.37 
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1 

‘lA’iL.1-  21  l0.. 

GUl'DJUJ, 

HEAT  COPttU'JTICt;  EFFECT 

i 

i 

,  Cr/C[aj[n 

=  10. 0  , 

(hA)*  3  As  * 

=  1.0 

e  70 

Conduction  effect ,  —  7o 

fJTUo 

s 

o 

0.01 

0.02 

0.04 

0.03 

o :  c 

0.32 

l 

56.42 

0.28 

0.53 

0.98 

1.74 

2.90 

4.37 

2 

i 

77.34 

0.37 

0.71 

1.35 

2.49 

4.36 

7.06 

3 

87.30 

0.38 

0.74 

1.43 

2.68 

4.82 

8.03 

4 

92.60 

0.35 

0.70 

1.36 

2.60 

4.80 

8.29 

3 

95.59 

0.21 

0.62 

1  .22 

2.39 

4 . 53 

8.1 

10 

99 . 62 

0.10 

0.22 

.49 

1.14 

2.68 

5.96 

20 

99.99 

0.01 

0.01 

0.05 

0.23 

1.02 

3.65 
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*  r  4  fy  j  ij  ^ 


1 

\ 


C  fr%  =  r»  «; 

win'  t.iax  U*J 

Cr  /Crain 

-  1.0  , 

(hA)  -  0* 

-0.5 

« % 

Conduction  effect ,  ^r— % 

IJ  i' 

r- 

_ 

o.oi 

0.02 

0.04 

O.Q3 

0.16 

0.32 

1 

50.5? 

!  i 

0.2C 

i 

53 

0.°rf 

1 .  ?’ 

**  t  c 

fa  »  >  J 

3.99 

2 

65.65 

0.35 

.69 

1.32 

2.40 

6.1  3 

6.50 

*> 

7  2. 44 

1  l 

0.30 

1 

.78 

l  .49 

2.78 

4 . 92 

8.04 

4  i 

76.43 

0.*2 

.87 

1.6? 

3.04 

5.46 

9.04 

5 

79.10 

1 

! 

0.45 

! 

0.89 

1.72 

3.  ?e. 

5.86 

9.79 

10 

! 

85.44 

i 

0.56 

1.09 

2.12 

4.00 

7  8 

11.95 

20 

89.72 

i 

0.74 

1.44 

2.7° 

4.95 

8.58 

13.90 
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i 
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1  j  # 

-> 

9  ^  O  i  1 1  » ] 

-  0  >  0 

i*.  .  .'.  1 

V.  1  Fl  ^  *'  ‘C* ' 

\ 

Ui. 

»  rlu 

>  V 

**  ••  i 

(U) 

'  0,5 

-  % 

To  ’duct 

ioa  fj  tft 

Jo 

--  ■■  -., 

- 

{■  * 

*  K* 

o,C2 

0,c*4 

0 . 03 

C  .lt> 

0.32 

1 

j 

— 

\t  J  l 

:  1 

I 

0.61 

1.71 

’  .!'*» 

<"».  -r  f 

•S .  i  •* 

4.54 

%) 

77. 72  | 

j 

04t.r 

C.33 

i  .:7 

I .  •’  J 

f  :  4. 

7.76 

, 

{*■*  27 

j 

0.37 

l .  6G 

3.1' 

4.61 

9.11 

1 

4 

37.57 

,  ] 

0.42 

0..2 

1.  52 

3.10 

5.71 

9.73 

p. 

^  ■ 

95.56 

1 

I  i 

0,37 

1 

[ 

3.74 

1.47 

2.90 

5.53 

o  7#. 

i 

10 

j  i 

99.61 

1 

0,1 

0.22 

0.43 

l .  54 

d.l  3 

^  0 

9^,99 

°*- 

! 

'!!  .  02 

«.0£ 

j.  43 

1.9? 

4.06 
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TABLE  24  LONGITUDINAL  HEAT  CONDUCTION  EFFECT 


=  0.5 

**-  A 

»  *  r/^Min 

-  1.0  , 

(hA)  =  AS  ” 

0.25 

1 

«  7o 

Conduction  effect  % 

NTUoV  ' 

\! 

0 

.01 

.02 

.04 

.on 

.16 

.32 

1 

j 

1 

50.21 

0.34 

C .  65 

1.13 

2.01 

3.12 

4.!6 

2  ! 

64.39 

0.30 

0.75 

1.42 

2.61 

4.47 

6.89 

3 

70.90 

0.39 

0.77 

1.51 

2.89 

5.19 

8.38 

a 

74.91 

0.39 

0.79 

1.58 

3.09 

5.70 

9.43 

5 

77.57 

0.40 

0.81 

1.64 

3.26 

6.11 

10.24 

i 

10 

04.12 

0.49 

0.99 

1.99 

3.97 

7.53 

12.71 

20 

88.72 

0.66 

1.30 

2.54 

4.92 

9.07 

14.97 
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TABLE  25  LONGITUDINAL  HEAT  CONDUCTION  EFFECT 


t  Cr^Omin 

=  10.0  , 

*  * 
(hA)  =  An 

«  0*25 

e  7c 

Conduction  effect^  % 

HTUoN^ 

0 

•  01 

.02 

.04 

.08 

.  •  6 

.32 

1 

56.40 

0.43 

0.81 

1.47 

2,51 

3.92 

5.46 

2 

77.30 

0.58 

1*11 

2.10 

3.79 

6.36 

9.59 

3 

87.24 

0.60 

1.18 

2.27 

4,24 

7.42 

11.66 

4 

92.53 

0.57 

1.12 

2.22 

4.28 

7.77 

12.62 

5 

95.53 

0.50 

1.02 

2.05 

4,10 

7.73 

12.99 

10 

99.60 

0.18 

0.41 

1.01 

2.59 

6.18 

12.19 

20 

99.99 

0.01 

0.04 

0.22 

1.13 

4.29 

10.59 
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NTUo 


APPENDIX  I. 


Derivation  of  iterative  equations. - 

Equations  (1)  ami  (2)  from  page  c  can  be  put  in  the  following  form: 

t  —  f  -Cl  .  T  —  T  — T  (l.l 

*(•«  x's-»  +  j  1  '-'x’.jj  'H'-*1  .})  *v.<.  j;  V.  j.') 


"  v<-ti  y, 


V-  ^  ) 


(1.2) 


where 


N  I 


J 


I 


/  \ 


and 

V— * 

Equations  (1.1)  and  (1.2)  may  be  written  as 


(1.3) 


and 


r  U 


(1.4) 
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Dy  skiing  to  equation  (1.3)  the  product  of  (l-rC^,  )  (l  .4)  ana 


-i 


(*> 


vhere 


and  - 

-  /-* 


From  equation  (1.1) 


(6) 


where 


ft 


and 


For  the  elements  from  the  first  and  last  row  equation  (1.2)  becomes 


and 


(1.5) 


(1.6) 

\  .  -  i 
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respectively.  Tw  t  *e  same  procedure  ns  before  mid  solving  for  the  matrix 
outlet  temperatures,  the  following  equations  are  obtained: 


-4 


(3) 


and 


N 


(5) 


for  the  first  and  last  row  eLements  respectively,  where 


and  \J  - 

<r 


Equations  (7)  and  (8)  ntay  be  written  as 


(1.7) 


and 


(1.8) 


Equation  (1.8)  reduces  to 


52 


t  -i 
^w-,S)  Wi,a) 

and 


U  ‘3  /  'd' 


+  t3 


’  '-n  _‘y^  _-p 

(^’  ;})  ^  U^r,3*H)  °i^t*“,5)  ^Wr“  ,g+l)j  (1.9) 


-  I  , 


(1.10) 


for  the  first  and  last  row  elements  respectively,  where 


-  N  I 


-  V 


By  comparison  of  equations  (1.7)  and  (1.8)  with  equations  (1.1)  and 
(1.2)  it  follows  that  the  form  is  the  same  and  only  the  coefficients  are 
different.  The  same  holds  true  for  equations  (1.9)  and  (1.10)  with 
equations  (1.6)  and  (1.5)  respectively.  By  substitution  the  following 
equations  are  obtained: 


o  .  ..  V; 


’  <  h  itf.  i) 


(9) 


'.Nr’’ 


i  : 


I.  'I. 


(10) 


I  • 1 


,  1 


(U) 


for  the  first,  middle  and  last  row  elements  respectively  and 
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(12) 


which  is  the  same  for  oil  the  elements,  where 


I 


i 


; 


and 


54 


the  relieving  form: 


Convt: egsucc  cot*  Ji  t  ion ,  - 

Lquation  (i  *2)  car  b~  pat  in 


(l.n) 

n  the  energy  balance  equations,  the  aide  ot  Crn3x  was  assumed  to  be  the 
hot  side  so  that  by  the  second  law  of  thermodynamics 


and  from  equation  (1.11) 

4  i  .  .  — 


Rear?.*,  tv  i  *r  equation  (1.12) 


(112) 


(1.13) 

Again  by  the  second  law  oi  thermodynamics  the  matrix  temperature  variation 
in  the  direction  of  fluid  flow  is  ouch  that 


(l.U) 

<(  ^ 
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aUd 


(1.15) 


so  that 


(1.16) 


Substituting  equation  (1.16)  in  equation  (1.13) 


(1.17) 


Three  cases  are  investigated: 

1. -  If  C3  C,  then  equation  (1.17)  becomes 

>  — 

j  * 

2. -  If  Cj  »  0  then  equation  (1.17)  becomes 

y  j  ■  • 

3. -  If  C,  •  2C^  *  0  then  equation  (1.17)  becomes 


From  the  above  cases  it  can  be  seen  that  a  sufficient  condition  for 
equation  (2)  to  converge  is  that 


Substituting  the  values  for  the  constants  and  rearranging  yields 
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For  the  first  and  last  row  elements  the  convergence  criterion 


becomes 


which  for  substitution  of  the  parameters  gives 


This  condition  is  less  stringent  than  for  the  middle  row  elements. 

Similarly,  Cor  the  side  o.  Crain  the  sufficient  condition 
for  the  middle  row  elements  becomes 


which  for  substitution  of  the  parameters  gives 


For  the  first  and  last  row  elements  it  is  found  that 


which  requires  that 


is  the  sufficient  condition.  This  again  is  a  less  stringent  condition 
than  that  imposed  by  the  middle  row  elements. 
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r» i i  (even t  ial  equ.it  io  u .  - 


Hfius*  n  4  hi*s  given  a  theory  for  regenerator*  without  heat 
conduction.  When  the  wiit  ussumpt ivr.tt  -re  nr.  Se,  but  tin.  longitudinal 
conduction  of  heat  is  included,  (See  page  *  )  the  differential 
equations  obtained  are  _  : 

for  the  aide  of  Cmax  (gas  flowing  from  ;  -  0  to  *  L) ; 

i 


and  for  the  side  of  'reversed  gas  flow) 


A  temperature  scale  can  be  used  for  which  the  gas  entrance 
temperature  at  the  side  ol  CT1jn  is  zero  and  that  at  the  side  of  C 

MS 

is  unity.  The  boundary  coi.a  tions  are  then 


and 


APPENDIX  2 


Computer  program. - 

The  Fortran  system  was  used  to  program  the  problem  for  the  CDC  IbO'f 
digital  computer. 

The  sij:  dimensionless  parameters  are  used  as  input  together  with 
increments  and  factors  to  change  the  values  of  NTUo,  and  As,  so  that 
solutions  can  be  generated  for  several  combinations  of  the  parameters 
without  having  to  feed  into  the  computer  these  different  sets  of  para¬ 
meters  every  time. 

Another  input  is  the  number  of  subdivisions  of  the  three  streams 
together  with  the  increments,  so  that  when  the  computer  finishes  the 
first  series  of  runs  for  a  given  number  of  subdivisions  it  automatically 
increases  the  subdivisions  by  £he  specified  increments  and  the  two  values 
of  effectiveness  obtained  are  extrapolated  to  an  infinite  number  of 
subdivisions. 

In  order  to  obtain  the  initial  estimates  of  matrix  temperatures, 
the  problem  is  first  solved  neglecting  the  longitudinal  heat  conduction 
as  explained  previously.  For  this  case  the  equations  derived  Ln  Appendix  1 
are  simplified  to  (for  detailed  derivation  see  Reference  1) 
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The  inlet  matrix  temperatures  necessary  to  initiate  the 

problem  for  zero  heat  conduction  are  estimated  by  assuming  \,  , 

(See  1'ig.l'.  )  rnd  feeding  this  value  into  the  computer.  The  retraining 
inlet  matrix  temperatures  are  obtaLned  from  the  relation 

H‘  ij  ~  'i  ' 

The  number  of  pusses  necessary  to  meet  the  desired  accuracy  in  the 
solution  changes  accordingly  with  the  initial  estimates  of  matrix  inlet 
temperatures,  but  considering  that  the  computer  time  required  to  work  the 
problem  for  zero  heat  conduction  is  so  small  as  compared  to  the  case  with 
conduction,  that  these  estimates  do  not  become  important.  Only  for  high 
values  of  Ci/Qain  these  estimates  may  be  of  significance  since  for  this 
case  the  change  in  the  reversal  condition  after  every  pass  is  very  small. 

In  general,  it  was  found  that  an  estimate'  T,  jv 


of 

0.4  to  0.6 

for 

^r  /^in 

<5.0 

and 

0.7  to  0.9 

for 

Cr/Cmin 

^  5.0 

together  with  the  above  relation  for  Tx,^  should  converge  to  the 
solution  in  a  reasonable  time. 

The  temperature  distribution  is  calculated  by  repetitive  use  of  the 
iteration  equations. 

The  effectiveness  may  be  defined  as 


0  Hi  ax  (  ^  XoJ  _  ~trio  ~ 

C-n viw  (t^xi  —  T-n,*)  Txi- Vi 
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irjio  '  s-  t  j  oiv;  of  i> 


-rcMcu 


therefore,  the  expression  for  effectiveness  for  the  definition  reduces  to 


N  , 


which  is  the  expression  used  for  the  computation  of  the  effectiveness. 
The  heat  balance  error  is  computed  from  the  relation 

(  A ,  -  u  ) 


error  = 


J.  V 


which  reduces  to 


error  = 
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The  extrapolated  of Cectiveneas  is  conput«d  fren  the  relation 


(iv 


2) 


i 


where 


The  input  parameters  tn  Fortran  notation  are 
^  =  Sain^aax  >  i  2  -  Cy/C,,^  ,  F3  =  (liA)  , 

PA  =  NTUo  ,  P5  «  \  ,  ?6  =  As* 

The  initial  values  of  NTUo,  \  and  As*  are  changed  by  the  following 
relations 

?A  «  p4  PAIN 
P5  =  P5'*?5FA  +  P5IN 
P6  =  P6  -  p6lN 

where  P4LN,  P5FA,  P5IH,  P6IN  are  also  input  values.  The  inputs  pAFI, 

P5FT.  and  P6FI  are  the  final  values  of  P4,  P5,  and  F6  respectively. 

The  program  is  written  so  that  the  computer  after  working  the 
initial  set  of  parameters,  i'  reduces  '6  by  P6IN  every  time  until 
P6  becomes  less  than  P6FI.  Then  it  sets  F6  equal  to  the  initial  value 
and  changes  P5.  The  procedure  is  again  repeated  until  P5  becomes  greater 
than  P5FI.  Then  P5  is  set  equal  to  the  initial  value  and  P4  is  increased 
by  P4IN.  The  cycle  is  repeated  until  P4  becomes  equal  to  or  greater  than 
F-4FI. 
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L/iiri r  f  ■  t  ,  a  •'•iz  ;  ^  trr^r  bv  t be  *  s  ^ t  Vr  z^ 

Ftr.fc  >  5v*.nt*i  vis*'  S  'ir^v  ':  to  wake  s !  *c  f.  U  v?rt«  »n  *J  v  Z  it:  *.  *“» 

r'  '■  or  d;yr,  no:  ^oiv  rf  wi.ro 'n  th*  i-inxiarjii’  iu,  v  •  -  « r»  ■  -n 

*)f»s  fii‘  ,;.  » ■$  *  i .  ?r.  7b  <?ho*>  the  co^ipi -1  e  fl c<f  f  v,  ?  f»- 

ore;;  »  v*hlch  l*  follow*!  by  the  For  tun  listing  of  the  t 


H  X.  J 
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I  1  INPUT  VALUES J 

M  =  PASS  NUMBER 

2  COMPUTE  INITIAL  Tx(j,l)  |  IC  =  1,  INITIAL  SUBDIVISIONS 

IC  S  2,  INCrtE  >;NTLD  SUBDIVISIONS 


[3  SET  =  1 

.0  AND  tn(l,g)  =  0  | 

NO  CONDUCTION  CASE 

COMPUTE  Tx(t,j)  ,  Tn(r,g)  ,e  AND  ERROR 

T 


’  44  SET  Afi"  *  INITIAL  VALUE 
INCREASE  X^. 


46  SET  \^=  INITIAL  VALUE 


FIG.  ?■ 


GENES AL  FLO*  DIAGRAM 


I  4  COMPUTE  r  '  NST  *N  IS 
'5  SET  K  c  1 

1  6  ITZ'Wi.V  'i  s.->s  or  C3flx 
7  ITERATION  ON  ST^S  CF  C^in 
J  6  COMPUTE  e  AND  ERROR 


ERROR  £  Si~ClFI  ED  ^XIMUM©> 


NC 


i  15  STORE  Tx(i , j ) »  Tn(f,g) 


©  f  1 7  SENSE  SWITCH  [  © - j 

[©  PRINT  e  AND  ERROR] 


CONDITION 
REPLACE  Tx(itl) 

Bif  Tn(f,Nn*l) 


<2> 


] — © 


22  INCREASE  Nr  ,  Nx,  Nn 


ILO'A  DI A  ».<AM  FOR  NO  CONDUCTION 
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FIG.  7b 


PROORAM 

remarks 

PROGRAM  BAHNKF 

01  MENS  JON  TH(4Q,40),UH(40,40I  ,10140,401  ,UC  (40,40)  ,  ERR  (92)  ,IC0l92>* 

1 UHC ( 40*401. UCCC 40, 4Q 1 , UH 1(401 ,UH2 (40),fF(20,20),FR'20,201 

)  FORMAT  1 6f 5.2 1 

RE  AO  l,Pl,P2,P3,P4,P41N,P4f 1 

2  FORM AT ( 6 1 5 1 

RE  AO  2»NM,NH,NC,NM1 , NHl,  NCI 

INPUT 

NM  *  Kr,  NH  -  Kz,  NC  -  Na 

S  FORMAT ( F  5 • 3 1 

READ  3, UHC 1,11 

E  M*NM 

UO  401  K®  1 » NM 

F=K 

401  UHCK, 11  =  1 (EM  +  l.O-E l•UH{l ,  1  )  1/EM 

INPUTt  0H(1,1)  -  Tx(l , I  ) 

ESTIMATES  OF  T*(ifi) 

4  FORMAT ( 7F  7 .2 1 

RE  AO  4, P5,P5IN,P5FA,P5f 1  , P6 , P6 1 N , P6f  l 

INPUT 

S5*P5 

S6  =  P6 

STORE  INITIAL  P5  AND  P6 

5  NM2=NM ♦ NM I 

NH2®NH*NH 1 

NC2®NC*NCI 

NUMBER  OF  SUBDIVISIONS  FOR  2nd  RON 

NE  1  -NM • ( NH  *NC 1 

RE  UNE  ( 

NE2=NM2» (NH2+NC2 ) 

RE2=NE2 

R£-REW1RE2~RE  1  1 

RE  •  VARIABLE  USED  TO  EXTRAPOLATE  EfTECTI VENESi 

EM2=NM2 

1)06  L=  1.NM2 

t  -t 

RT*1.0-t(E— 1.01/EM2) 

6  UH2C  L  1  =UH  (  I,  Il'KT 

ESTIMATES  Or  Tx(lfl)  TOR  2nd  RUN 

7  FORMAT!// I4HN0  OF  E L E MEN l S3X , 3 l 5 , 4 X , 3 ( 5 ) 

WRITE  OUTPUT  TAPE  2, 7  »NM,NH»NC,NM2, NH2  *  NC2 

PRINT  7,NM,NH,NC,NM2*NH2,NC2 

8  I  C=  1 

INDICATES  lit  RUN,  USING  INITIAL  No  Or  3UBDIV. 

9  EM=NM 

EH  =  NFf 

6C«NC 

riXED  TO  FLOATINO  POINT  CONVERSION 

NMM=NM- ( 

NHP’NHF ( 

NCP=NC+I 

NM  P  =  NM ♦  I 

t 

DEFINITION  OF  NEW  VARIABLES 

00  10  J* 1 ,NH 

10  TH (  1  *  J  1  =  1 .0 

UO  11  L-1.NC 

1)  T C  C 1 , L 1 -0. 0  .  .  

SET  >  1.0  AND  tn(l,g)  *  0.0 

) 2  A  1 0  = (.0*EM/(PI«P2»EH1*(2.0-EM*P31/(P4#PI«( 1 . 0*  P  3 1 ) 

A  1  =  2. 0/A  10 

A20= 1.0*1 p 1 • P2 • E  H 1 /EM+ (2 ,0«EH»P3*P2 l/(P4.(|.0*P31l 

A2=2.C/A20 

A30=1.Q*EM/(P2»FCI  ♦ ( 2 .0 • E M ) / l P4 • < 1 . 0 ♦ P3 1  1 

A3-2.0/A30 

A40= 1 .0*(P2»EC 1/EM  +(2.0*EC«P21/(P4»(  1  • 0*  P3 1 ) 

A4=2 . 0/ A40 

CONSTANTS  FOR  THE  NO  CONDUCTION  CASE 

C  HOT  SIOE  NO  CONDUCTION 

121  M  a  1 

M  DENOTES  PASS  NUMBER 

13  00  ) 5  J=I.NH 

DO  IS  1  =  1*  NM 

THClF(,Jl  =  THCItJUAI-CUH(l,Jl-IH(|,Jll 

15  UH ( I  *  J  *  1 1 =  UH ( I , J 1-A2»  CUHCI ,J|-TH(I , Jl ) 

ITERATION  ON  SIDE  Or  CB*x 

C  COLO  SIDE  NO  CONOUCTION 

00  17  K= 1  *  NM 

17  UC(K ,  1  1=UH(NM+  1-K*NH*  11 

UO  19  L  = 1 »  NC 

0019  K=(,NM 

TCCK*  1,L)  =  TC(K,L  1*A3»  CUC  (K.U-TC  (K.L  1  ) 

19  UC ( K , L ♦ 1 1 =UC ( K , L I-A4MUC 1K,L)-TC(K,1>) 

ITERATION  ON  SIDE  OF  C»iD 

01=0.0 

00  20  J= 1  *  NH 

20  GI>TH(NM«1.J1«GI 

THO=Gl/tH 

TH0  •  [tz(Krd,J)j  Brg 

G2=0 .0 
00  2  1  L  =  1 1 NC 
G2*TC( NM* 1 ,L I +  G2 
TC0lM)=G2/EC 


TCC  -  EFFECTIVENESS 


ERR(M1=1.0-IP1«TCO!MI )/( 1.0  Til  I 
IFIIRR(M)  )23,24,?4 


COMPUTE  HEAT  BALANCE  ERROR  AND  COMPARE  IT 


c««ini*-cnn»ni 

24  IE(ERR(M)-4*0f-6l  30, 30 ,25 

WITH  MAXIMUM  SPECIFIED 

1  2S  (FI90-MI  26,26,215 

j  26  IF  1  ERR! M- 1 1-ERRI M 1)  2 9,29,2? 

|  27  PAUSE  1 1 

CHECK  CONVERGENCE 

1  28  IFI SENSE  SWITCH  1 1  32, 121 

PAUSE  11  =  SOLUTION  CONVERGES 

29  PAUSE  12 

PAUSE  12  =  SOLUTION  DIVERGES 

GO  TO  28 

215  DO  22  1=1, NM 

22  UH 1 I ,  1 1 =UC 1NM+  1-ItNC+l  ) 

REVERSAL  CONDITION 

M-Ht  1 

GO  TO  13 

301 

302 


T 2  =  TCO l M- I )- ICO ( M  J 
1FIT2)  301,302,302 
T2= - I 2 

IF  (  T2-4.0E-61  305,  305, 2 ‘j 


305 

31 

312 


1  F  I  (C-  I  1  31*31,  33 
DO  312  J=1.NM 
UH  1  I J  )=UH( J ,  1  1 
6  F  C  1  ,  1  1  =  TCOl M 1 
ER  *  1,  I  )  =  ERRt  M 1 
GO  TO  3V  _ 


CHECK  DIFFERENCE  IN  EFFECTIVENESS  IN 
THE  LAST  TWC  PASSES 


STORE  Tx(itl)  ,  EFFECTIVENESS  AND 
ERROR  FOR  THE  1st  RUN 


FORMAT ( 4FB.2,  1 5, F 10.6, F9.2, »  31 

PRINT  32,PI,P2,P3,P4,M,TCO(*>,ERklM) ,!C 

P^USE  13 

IFCSENSE  SWITCH  1)  32 1 , 1 2 1 
IF  1  SENSE  SWITCH  2)  1,322 

NM=NM2 

NH=NH2 

NC=NC2 

G0_T0  _5.. _ _ 

EFF=TC0|M1* 

B 1  =  t FF 


INCREASE  NO.  or  SUBDIVISIONS  AND  REPEAT 
1st  RUN 


l.O+RE  1-EFI  I  ,  I  l«Kh 


1  34 

FORMAT! //4F8. 2, 

WR1IE  OUTPUT  T l 

00  35  L= 1 , NM 

35 

UH2 I L 1 =UH1 L , 1 ) 

39 

I  A=2 

I  B=  1 

+  9. 21  ,F  10 .6) 

/  ,P3,  P4  ,EF  I  1  i 


EXTRAPOLATE  EFFECTIVENESS  AND  STORE 


1)  ,  E  K  (  1,1)  ,  Tl  0!*>  ,  E*R(*>  , 


STORE  Tx(iti)  FROM  2nd  RUN 


40  DO  44  1=1, NM 

Ml  00  42  J-  1.NHP 

42  UHC ( I , J 1=UH( I , J ) 

43  00  44  1=1, NCP 

44  UCC!  1 ,1  )=UCI  1,L  1 

C  CONSTANTS 

50  £1=(P4*1  1.0+P3) 1/12. 0*EM  ) 

61 1* I.0+E1 

F  2=  C P2«EC)/EM 

E  3=0 . S*PS*EM 

E4=E 1+{ E 1 1*1 2.0*E3+E21 ) 

E5  =  E  HIE  1  I •( E2  +  E  3)  1 
F 1=1 2.0*6 1I/F4 

F2=( E 1  +  1  El  1*1 2.0*E3-E2 ))  1/EM 

F  3* 1 E  3*  E  1  11/EM 

F4  =  1  1.0-E I  1 / 1  1  1 

F5  =  E  I/El  1 

F6=( 2.0*E  1  1 /t 5 

F7=1E  l+l  E  I  1*(E3-E2  )  1  1/kS 

F8=  (  E  3*6  1H/P5 

C1=1P1*P4*1 I.0+P3I J/l?.0*P3*tM) 

Cll=1.0+C1 

C2=1P1*P2«FH)/EM 

51  C3-1P1*PS*EM1/I2.0*P61 
C4*Cl  +  (CU*(C2  +  2.0*f  3)  > 

CS=C  !♦( C  1 1 « I C2+C  3 1  I 


SUBSCRIPTS  FOR  ERROR  AND  EFFECTIVENESS 


UHC(I.J)  =  Tx(i|J>  ,  UCC(I,L)  .  Tq{ f f g ) 
FOR  CONDUCTION  CASE 


CONSTANTS  FOR  CONDUCTICN  CASE 


5  3 
C 

60 

6) 

62 


64 

66 


67 

60 

69 
C 

70 

71 

72 


73 

74 


75 


D1*(2. 0*C1  )/C4 

D2=IC1*(C1  Nt2.0*C3-C2)  )  J/C4 

03* I C  3*C  1  I  )/C4 

04= (  I  . 0-C  ll/CII 

05  =  C  1  /C  I  1 

06= I  2 «0*C  I  1/C5 

07= I C  I  ♦ I C  1 1 *1 C3-C2 1  )J/C5 

06= ( C  3*C  1  I  1/C5 

M=  I 

HOT  SIDE  WITH  CONDUC  1  I  0»M 
00  67  J»1,NH 
L=  ) 

UHCt  l ,  J+  1  )  =1)6  •  T  H I  1 ,  J  )-0  7*UHC  I  1  ,  J  )  *00  *  (  UHC  I  ?  ,  J  1  ♦  UHC  (  2  ,  J  *  I  1  » 

TH( 2, J )=04*THt  I,  J  )  *05*  (  UHC  (  I  ,  J  1  *UHCI  I  ,  J*  Ml 
0064  1=2, NMM 

UHC  (  1  ,  J*1)=U)*THt 1  ,  JI-02HJHCI  1  ,  J  I  *03  •  (  UHC  (  1-  1  ,  J  )  *UHC(  )-) , J* 1  1* 
l  UHC  (  1  *  1  , J 1 *UHC I  IH.Jtlll 
TH(  ) *  I , J 1=04*TH( I , J ) *05* l UHC t I  , JJ+UHCI I , J* III 

UHCINM,  J*  1  )=06*TH(  NM,  J  )-C7*LHCINM,J)  ♦ D8 • I UHC l NMM , J  )  *UHC(  N*»m,  J*  1  )  > 
THl  NMP,  J»*04*THlNM,  J  J  +D5MUHCINM,  J  )*UHCINM,  J* 1  J  ) 

L  =  l*  1 

) F ( 3-L )  67,67,61 

CONTINUE  _ _ _ 

DO  69  K=),NM 

UCC1K,  I  1 =UHC ( NMP-K , NHP ) 

COLO  SIDE  WITH  CONDUCTION 

00  75  L=1,NC 

J=1 

UCCI  1,1*1  1 =F  6*TC I  1  ,  L  )  -  F  7  »UCC  1  1  ,U  *FH  •  (  UCC  I  2  ,  L  )  *UCC(2  ,  L  +  1  )  1 
TCI  2,l)=F4*TC<  1,  L1+F5MUCCI  1,L)*UCCI  1,1*111 
00  73  K  =  2 »  NMM 

UCCI K ,  L  ♦  1  I =F 1*TC(K,L)-F2*  UCC I K , L 1 ♦ F  3  *  I UCC I K- I ,L>  *UCC(K-I ,L* 1  )* 
1UCC IK*  I , L 1 *UCC I K*  I  ,L  +  I  1  1 
TCIK*  1 , L  )  =  F4*TCI K, L 1  * F 5* I UCC I K , L 1 ♦ UC C I K  ,  L *  I ) 1 

UCCI NM,L* l J  =  F6*TCI NM,L )-F  7*UCC I NM , L ) *f  8* I UCC I NMM, L 1 *UCCC  NMM, i  *1)1 
TCtNMP.L }=F4»TCI NM , L ) + F5 • I UCC I NM , L 1 + UCC I NM , l ♦ I) 1 
J=J*  I 

I  F ( 3- J 1  75,75,71 
CONTINUE 
G  1  =0  «  0 


DO  76  J= 1 , NH 
76  G1=THI NMP, J 1*G 1 
TH0  =  G1/EH 

G2  3  0  •  0 

00  70  L  = 1 , NC 


70  G2=TCINMP,L 1 *G2 
TCP  I M 1 =G2/FC 

ERR(M|= 1.0-1 P1*TC0( Ml )/(  I.O-TllO) 
1 F | ERR  I M ) )  00,00,01 
00  ERR1M)=-ERRIM) 

61  | F I  ERR  I M )-4 .  OE-6 )  90,90,  8? 

02  IFI90-M)  03,03,05 
83  IFIERRIM- 1 )-ERRtM» J  04,64,87 
04  PAUSE  22 

04  1  I F I  SENSE  SWITCH  II  042,53 

042  1FISENSE  SWITCH  2)  00,084 

85  M=H*  I 


00  06  1=1, NM 

06  UHC I  I  ,  I  1 =UCC l NMP- I , NCP  I 

_ GO  TO  60 _ _ _ _ 

07  PAUSfc  21 
GO  TO  04  I 

80  1  F  I  IC-  11  001,001,082 

08  I  EF|  I  A ,  181=0,0 
ER(  I  A,  101=0.0 
GO  TO  97  _ 

002  FORMAT |  1 1HS0L  OIV  FOR  I  3 ,  I 0x  ,  E 9 . 2  ,  F 0.  2 , 2  I  F 1 0 . 6 , F 9 . ? 1 1 

WR1TF  OUTPUT  1  APE  2 , 882  ,  IC  ,  05,  P6  ,  F  FII  A  ,  I  8  J  ,  PR  (  1  A  ,  l  f!  1  ,  TCO  I  M  )  , 
I _ QO  TO.  97 _ _ _ 


t  RR(  Ml; 


CONSTANTS  FOR  CONDUCTION  CASE 

_  J 


iteration  on  SIDE  OF  C..X 
TWO  PASSES  PER  COLUMN 


ITERATION  ON  SIDE  Or  Cain 
TWO  PASSES  PER  COLUMN 


884  FORMAT  I 4F 8 ♦ 2, E 8. 2, F 8.2,2*, |$VFI0.6,E9.?1 

PRINT  884,  P  I,P2.p3,P4,PS,P6, 1C.TC0IM  ]  ,  F  RR  (  M  1 
PAUSE  23 

IFISENSE  SWITCH  1)  BPS. S3 

BBS  IFISENSE  SWITCH  2)  114, 8e 

89  IFIIC-11  891,891,93 

8  9  I  EFUA.  I  8  I  =  TC  0 1  «  1 

ER I |A, I B 1 =ERR ( M 1 

GO  TO  102  _ 

90  TI*TCOIM-l 1-TC0IMI 
I F | T  II  9  1 ,92,92 


T  1=-  T  I 

IF ! T 1-4 • OE-6 I  89,89,82 


EFF  =  TCO(MI*(  I.O*REI-FFIH,IM»Rb 
82*  C I B 1-6FF ) • 100.01/P 1 


97 

971 

972 

973 

974 


FORMAT  1 32 X, E8.2.FB.2.2IF I0.6,E9.?I,F10.6,F8.31 

WRITE  OUTPUT  TAPE  2, 94 , PS , P6 , E f 11  A , 1 8 1 , fc R I  I  A , I  8 1 , TCO I M ) , E RR l M I , fc FF 
1,82 

IF  I  SENSE  SWITCH  3)  9S,  102 

FORMAT  I F 8 . 2 , E9.2.FB.4 1 

PRINT  95 , P4, P5 , 82 

GO  TO  102 

00  974  I=),NM 

00  972  J= I , NHP 

UHC ( 1 , J )=UH{ I , J ) 

DO  974  L-  I  , NCP 
UCCt  1  ,U=UC(  I  ,  L  ) 


102  IFISENSE  SWITCH  3)  131,132 

131  FORMAT! /4F8.2,E8.2,F8.2, 3x , I  3/ !  1 SF  7. 4 1 ) 

WRITE  OUTPUT  TAPE  3 ,  1  3  1  ,  P  1  ,  P  2  ,  P  3  ,  P4,  PS  ,  ,  I  C  , 

UUHCt  J,  1  1  , UCC  (  NMP— J  ,NCP)  ,J  =  1  ,NM  1 

132  P6=P6-P61N 


18=  IB*  1 

1 F { P6-P6F 1 )  I03.103.S1 

103  P6-  S6 
_ I  B=  I 

104  PS*PS*PSFA*P5IN 


1  A= 1 A ♦ ) 

1FCP5FI-PSI  10S,S0,S0 

105  P5-S5 

1FIIC-1)  106,106,108 

106  IC=2 


NMFNM2 

NH*NH2 

NC*NC2 

DO  107  K=1,NM 

107  UHIK,  I  )=UH2IK) 

_  GO  TO  9 

108  I F 1 P4F I— P4 1  113,113,109 

109  P4=P4*P4 IN 


NM-NM-NM 1 
NH*  NH-NH  1 

NC  =  NC-NC 1 
00  111  K-  I  » NM 
111  UH1 K, I I-UHIIK ) 

__  GO  TO  6  _ _ 

113  PAUSE  SO 

IFISENSE  SWITCH  11  114,1 

114  DO  1  IS  L *  I  , NM 

115  UH ( L ,  I  1 -UH2 1 L 1 
GO  TO  5 

ENO 

ENO 


ro 


UNCLASSIFIED 
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